Acceleration techniques for semiclassical Maxwell–Bloch systems: An application to discrete quantum dot ensembles

The solution to Maxwell–Bloch systems using an integral-equation-based framework has proven effective at capturing collective features of laser-driven and radiation-coupled quantum dots, such as light localization modifications Rabi oscillations. Importantly, it enables observation the dynamics each dot in large ensembles a rigorous, error-controlled, self-consistent way without resorting spatial averaging. Indeed, this approach demonstrated convergence containing up 104 interacting dots (Glosser et al., 2017). Scaling beyond tests limit computational horsepower, however, due O(NtNs2) scaling (where Nt Ns denote number temporal degrees freedom). In work, we present algorithm that reduces cost analysis O(NtNslog2Ns). While foundations rely on well-known particle–particle/particle–mesh adaptive integral methods, add refinements specific transient with multiple derivatives. Accordingly, offer numerical results validate accuracy, effectiveness utility analyzing dots.